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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC)  for  the  Directorate  of  Technology. 
The  results  of  the  experiments  were  obtained  by  ARO,  Inc.,  AEDC  Division  (a  Sverdrup 
Corporation  Company),  operating  contractor  for  the  AEDC,  AFSC,  Arnold  Air  Force 
Station,  Tennessee,  under  ARO  Project  Numbers  V32I-P1A,  E32P-POA,  and  P32M-01. 
Tire  Air  Force  project  managers  were  Capt.  Stanislaus  L.  Ludwig  (CF)  and  Mr.  Elton  R. 
Thompson.  The  data  analysis  was  completed  on  August  13,  1979,  and  the  manuscript 
was  submitted  for  publication  on  September  20,  1979. 
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1.0  INTRODUCTION 

1.1  PURPOSE  OF  THE  EXPERIMENTS 

A  common  feature  of  many  combustion  flow  fields  is  the  recirculation  zone 
established  by  turbulent  mixing  in  a  constant-area,  axisymmetric  duct  of  two  coaxial 
flows  (Ref.  1).  An  example  relevant  to  aerospace  applications  is  the  so-called  “sudden” 
expansion  or  dump  1  combustors  used  in  ramjet-rocket  propulsion  systems  that  are 
designed  on  the  principle  of  establishing  and  maintaining  combustion  in  regions  of 
recirculating  flow  within  the  combustor.  Analytical  techniques  must  be  developed  and 
evaluated  for  analyzing  and  predicting  the  turbulent  mixing  and  combustion  phenomena 
for  various  configurations.  However,  there  is  a  lack  of  experimental  flow-field  information 
for  such  configurations.  Therefore,  the  primary  purpose  of  these  experiments  was  to 
obtain  data  to  aid  the  development  of  ducted,  recirculating  flow  models.  In  addition,  the 
laser-Raman  diagnostic  technique  was  to  be  used  to  obtain  the  data  to  demonstrate  the 
utility  of  the  technique  in  hostile,  AEDC  test  cell  environments. 

1.2  OBJECTIVES 

The  primary  objective  was  to  use  laser-Raman  scattering  to  provide  a  flow-field  map 
of  rotational  temperature  for  two  combustor  configurations.  Because  of  the  high  gas 
density  (near  atmospheric)  in  the  flow  field,  the  static  (or  translational)  and  rotational 
temperatures  are  essentially  equal.  The  use  of  a  pulsed  laser  with  a  nominal  pulse  width 
of  25  nsec  made  the  measurements  essentially  instantaneous.  Therefore,  a  precise  average 
temperature  could  be  obtained  by  accumulating  a  large  number  of  these  instantaneous 
measurements  at  each  flow-field  position.  A  second  objective  was  to  map  species  number 
density  for  the  two  configurations.  ••  * 

1.3  DESCRIPTION  OF  THE  FLOW  FIELD 

Figure  1,  reproduced  from  Ref.  1,  shows  the  essential  features  of  the  flow  field.  For 
these  experiments,  the  higher-velocity,  primary  jet  stream  was  air,  and  the  lower-velocity, 
secondary  stream  was  hydrogen.  Whenever  the  entrainment,  or  pumping  capacity,  of  the 
primary  jet  is  greater  than  the  mass  flow  of  the  secondary  stream,  the  eddy  of  the 
recirculation  region  is  generated  (Fig.  1 ).  The  recirculation  eddy  provides  mass  flux  to 
balance  the  jet  entrainment  requirements.  Although  the  recirculation  eddy  can  exist  with 
or  without  secondary  flow,  the  recirculation  eddy  structure  disappears  if  the  secondary 
mass  flow  rate  is  increased  to  some  critical  value.  The  duct  flow  will  then  be 
characterized  by  positive  axial  velocities  throughout. 
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2.0  EXPERIMENTAL  SETUP  AND  OPERATION 


2.1  TEST  CELL 

The  experiments  were  conducted  in  the  AEDC  Engine  Test  Facility  Propulsion 
Research  Cell  (RIA-1),  detailed  in  Refs.  1  and  2  and  shown  as  modified  for  these 
experiments  in  Fig.  2.  The  0.6-cm  (l/4-in.)-wide  slot  permitted  viewing  of  the  laser  beam 
focal  volume  over  the  full  width  of  the  flow  field.  Observations  were  made  through  a 
22.9-cm  (9-in.)-diam  quartz  window.  Entrance  and  exit  ports  were  added  to  the  top  and 
bottom  of  the  slot  to  permit  the  high-energy  laser  beam  to  be  focused  within  the  flow  V 
field.  These  ports  were  also  made  of  quartz. 

The  test  chamber  was  a  3.5-m  (S-ft)-long  stainless-steel  duct  with  an  inside  diameter 
of  13.3  cm  (5.24-in.).  A  mechanically  driven,  axially  traversing  nozzle  assembly  was 
mounted  inside  the  duct,  and  an  O-ring  system  provided  the  seal  between  the  assembly 
and  the  duct.  The  nozzle  assembly  included  either  a  1.33-cm  (0. 524-in.)-  or  a  6.65-cm 
(2.62-in. )-diam  nozzle  that  produced  the  primary  air  jet  and  an  annular  assembly  that 
injected  hydrogen  into  the  duct,  the  walls  of  which  were  cooled  by  an  external  water 
spray.  To  stabilize  the  flow,  the  secondary  hydrogen  flow  passed  through  two  porous 
plates  and  a  screen  pack.  A  spark  plug  was  located  downstream  of  the  observation  point 
to  ignite  the  combustible  hydrogen-air  mixture. 

2.2  LASER  SYSTEM 

The  laser  system  was  mounted  on  the  top  tier  of  an  aluminum  table  and  directly 
above  the  spectrometer  system,  which  was  mounted  on  the  lower  part  of  the  table  (see 
Fig.  3).  The  aluminum  table  was  securely  clamped  to  a  hydraulically  operated  traversing 
table  that  was  used  to  move  the  entire  system  vertically  and  horizontally. 

The  laser  was  a  Holobeam  600  series  pulsed  ruby  operated  in  the  Q-switch  mode. 
The  ruby  rod  was  15.2-cm  (6  in.)  in  length  and  1.2  cm  (0.5-in.)  in  diameter.  A  Pockels 
cell  provided  the  Q-switching  mechanism;  the  pulse  repetition  rate  was  15  per  minute. 

The  output  beam  from  the  laser  was  expanded  in  the  ratio  3:1  and  focused  with  a 
1000-mm- focal-length,  fused  silica  lens  at  a  point  inside  the  test  chamber.  After  the  laser 
beam  passed  through  the  chamber,  a  lens  reduced  the  size  of  the  beam  diameter  to 
permit  detection  of  the  beam  energy  by  a  pyroelectric  detector.  Appropriate  beam 
attenuators  were  installed  at  the  entrance  to  the  pyroelectric  detector  to  prevent  damage. 

A  second  type  of  laser  energy  monitor  was  installed  at  the  rear  of  the  Pockels  cell. 
A  portion  of  the  very  weak  beam  transmitted  by  the  rear  mirror  of  the  Pockels  cell  was 
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split  off,  passed  through  laser  line  filters,  and  detected  by  a  photomultiplier  tube.  The 
output  of  the  photomultiplier  could  be  either  displayed  on  a  storage  scope  for  inspection 
of  laser  pulse  shape  or  integrated  to  give  a  measure  of  the  laser  output  energy. 

To  facilitate  alignment  of  the  ruby  laser  system,  a  helium-neon  (He-Ne)  laser  was 
installed  behind  the  rear  laser  energy  monitor.  The  He-Ne  laser  was  initially  aligned 
through  the  center  of  all  the  ruby  laser  optics  to  keep  the  center  of  the  ruby  beam 
coincident  with  that  of  the  He-Ne  beam.  Any  necessary  adjustments  or  periodic  checks 
on  system  alignment  could  be  accomplished  by  using  the  He-Ne  beam. 

2.3  SPECTROMETER  SYSTEM 

Two  25.4-cm  (1 0-in.)-focal-length  quartz  lenses  were  used  to  collect  the  laser 
scattered  radiation  and  to  image  the  observation  volume  with  unit  magnification  onto  the 
2-mm-wide  entrance  slit  of  the  spectrometer.  An  aperture  at  the  entrance  slit  selected  a 
0.3-cm  (1/8-in.)  beam  length  for  observation. 

The  spectrometer  was  an  0.85-m-focal-length  double  spectrometer  with  1200  g/mm 
gratings.  The  regular  single-exit-slit  assembly  was  removed  and  replaced  with  a 
2-mm-wide,  double-slit  assembly.  Exit  slit  radiation  was  coupled  to  dual  detectors  by 
means  of  adjustable  mirrors  and  lenses.  The  detectors  were  cooled  (-26° C)  RCA® 
C31034A  photomultiplier  tubes  (PMT). 

2.4  DATA  ACQUISITION  SYSTEM 

The  technique  used  for  this  application  integrated  over  time  the  PMT  output  current 
for  an  interval  no  less  than  that  of  the  pulse  duration.  For  background  radiation 
subtraction,  a  similar  integration  was  performed  without  laser  excitation.  Because  of  the 
very  short  pulse  duration,  extremely  fast  integration  was  required.  The  LeCroy  Model 
2250L  integrator  was  found  to  be  adequate  for  this  application.  Each  channel  had  a 
common  gate  control  and  individual  analog-to-digital  (A/D)  converters.  The  converter 
outputs  were  stored  in  a  32-deep,  first-in- first-out  (FIFO)  memory;  consequently,  up  to 
16  laser  shots  were  possible  for  a  single  run  series.  These  results  were  averaged  and 
printed  out  as  individual  laser  shots  and  averaged  for  the  run.  The  pyroelectric  detection 
system  was  also  read  by  the  processor  and  similarly  averaged. 

The  block  diagram  of  the  system  is  shown  in  Fig.  4.  The  corresponding  timing 
diagram  is  shown  in  Fig.  5.  The  data  interval  sequence  began  upon  generation  of  the 
O-switch  trigger  pulse  which  occurred  after  the  flash  lamp  excitation.  This  pulse  occurred 
610  Msec  after  the  flash  lamp  excitation,  which  corresponded  to  optimum  laser 
performance,  and  approximately  250  Msec  before  the  energy  pulse.  The  charging  of  the 
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laser  capacitor  was  controlled  from  the  processor  and  was  repeated  a  preselected  number 
of  times.  The  processor  counted  the  interval  flags  and  compared  the  count  to  a  preset 
number  in  memory.  The  laser  was  in  a  free-running  mode  and  fired  automatically.  The 
FIFO  memory  was  read  out  to  the  processor  at  the  end  of  each  series.  The  Q-switch 
synchronization  (sync)  pulse  triggered  a  storage  oscilloscope  that  was  used  to  monitor  the 
timing  features  of  the  data  acquisition  functions.  The  oscilloscope  gate  pulse,  in  turn, 
triggered  the  gate  pulse  generator.  System  delays  and  an  adjustable  delay  in  the  pulse 
generator  accounted  for  the  250-nsec  separation  between  the  sync  pulse  and  the  Q-switch 
energy  pulse  and  lor  any  additional  delay  differentials  in  the  interconnecting  coaxial 
lines.  The  adjustable  delays  were  selected  to  insure  coincidence  of  the  data  gate  and  the 
detector  signals  at  the  integrators.  Differentials  in  lengths  of  cable  and  light  paths 
between  the  laser  energy  monitor  and  PMT  signals  were  accounted  for  by  adjustments  in 
coaxial  line  length.  These  adjustments  were  made  by  applying  appropriate  signals  to  the 
storage  oscilloscope  and  adjusting  the  delay  on  the  gate  generator  while  making  repeated 
laser  firings. 

The  oscilloscope  gate  was  used  to  trigger  a  one-shot,  multivibrator  control  pulse 
(Figs.  4  and  5),  The  data  interval  pulse  (OS/1)  is  a  flag  used  by  the  computer  to  indicate 
that  data  input  to  the  integrator  has  been  started  and  completed.  Background  data  were 
obtained  by  the  double-pulse  mode  selection  on  the  gate  pulse  generator.  Pulse  separation 
was  adjusted  to  12  ,usec.  Since  9  ,usec  are  required  for  conversion  and  1  m sec  for  transfer 
to  FIFO  memory,  the  overall  effect  was  to  follow  the  laser  scattering  data  with  a 
background  measurement  of  the  same  integration  interval.  After  N  =  1 6  background  pulse 
and  data  combinations,  the  sequence  was  shut  down.  The  data  interval  (OS/1)  flagged  the 
processor  either  to  reset  the  laser  charge  condition  or  to  terminate  the  run,  read  the  data 
from  the  FIFO  memory  of  the  integrator  module,  converted  the  results  to  decimal,  and 
printed  the  results. 

3.0  LASER-RAMAN  MEASUREMENT  TECHNIQUE 
3.1  SPECTRA  AND  INTENSITY  RELATIONS 

As  a  result  of  inelastic  collisions  between  laser  beam  photons  and  gas  molecules, 
laser  radiation  is  scattered  from  the  beam  at  frequencies  other  than  that  of  the  laser.  This 
phenomenon,  known  as  Raman  scattering,  is  detailed  in  Refs.  3  through  6.  Radiation 
scattered  at  longer  wavelengths  (lower  frequencies)  is  commonly  referred  to  as  Stokes 
Raman  scattering;  only  the  Stokes  wavelengths  were  used  in  these  experiments.  If  the 
inelastic  collisions  between  the  laser  photons  and  gas  molecules  cause  only  rotational 
energy  level  changes  in  the  molecules,  pure  rotational  Raman  scattering  results  and  the 
spectra  will  be  observed  in  the  <l000-cnr]  region  from  the  laser  wavelength.  For  the 
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flow  field  under  study,  the  gas  species  expected  in  the  spectra  were  N2,  02,  H2,  and 
H20.  With  the  Raman  Spectral  Program  (RASP)*  variation  of  the  relative  intensity  of  the 
spectral  lines  of  N2,  02>  and  H2  with  temperature  can  be  readily  calculated.  Figures  6 
through  8  are  calculations  of  the  Stokes  rotational  Raman  spectra  of  N2  and  02  in  air, 
and  the  spectra  have  been  convolved  with  the  transmission  function  of  the  spectrometer 
used  in  these  experiments.  The  variation  of  the  spectral  profile  with  temperature  can  be 
observed  clearly  in  Figs.  6  through  8.  Similar  calculations  for  a  mixture  of  N2  and  H2 
([N2]/[H2]  =  0.667)  are  shown  in  Figs.  9  through  11.  The  temperature-dependent 
intensity  ratio  R(Tr)  for  the  spectral  locations  7010  A  and  7061.2  A  defined  by 


I  (7061.2  A) 

R(Tr)  =  -21 - —  (1) 

Ic(7010  A) 

has  been  computed  using  RASP  and  is  shown  in  Fig.  12  for  both  air  and  pure  N2  for 
300  K  <  TR  <  3000  K.  For  the  flow  field  under  study,  the  actual  ratio  must  lie 
somewhere  between  these  two  calculated  extremes.  It  has  been  assumed  that  the  weak 
and  widely  dispersed  pure  rotational  spectrum  of  H20  can  be  neglected;  the  spectral 
measurements  of  Ref.  7  justify  this  assumption.  The  intensity  ratios  of  the  J=1  (7238  A) 
and  J=2  (7359  A)  H2  lines  have  also  been  calculated  and  are  shown  in  Fig.  13.  By  using 
the  measured  values  of  the  intensity  ratios  and  the  predicted  curves  of  Figs.  12  and  13, 
temperature  values  can  be  determined  from  both  the  N2/02-  and  the  H2-rotational 
Raman  scattering  profiles. 

The  intensities  of  the  Raman  lines  of  the  molecules  are  directly  proportional  to  the 
number  densities  of  the  molecular  species.  The  intensity  of  the  J=1  (7238  A),  H2 
rotational  line  can  be  expressed  as 

im(j=i,  h2)  =  n(u2)  [cf(h2i  Cp  m]-1  (2) 

in  which  n(H2)  is  the  hydrogen  number  density,  CF(H2)  is  a  proportionality  factor,  and 
Cf3(T)  is  a  temperature-dependent  term.  A  RASP  calculation  of  Cf3(T)  is  given  in  Fig. 
14.  If  it  is  assumed  that  temperature  has  been  measured  and  that  CF(H2)  has  been 
determined  by  calibration,  then  n(H2)  can  be  determined  from 

n(H2)  =  Cf(N2)  Irln( )  =1,  H2)  (3) 


*  Developed  at  AEDC  by  \V.  D.  Williams,  et  al. 
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Similarly,  the  number  density  of  air,  n(air),  and  N2,  n(N2),  can  be  determined, 
respectively,  from  the  following  relations:  1 

n(air)  =  Cp(air)  Cp  (T)  U  (7010  A),  (4) 

n(N2)  =  CF(N2)Cp  (T)  U(7010  A)  (5) 


The  results  of  the  RASP  calculations  of  CFl(T)  and  Cp2(T)  are  graphed  in  Fig.  15. 

Since  the  inelastic  collisions  between  laser  beam  photons  and  molecules  result  in 
vibrational  level  transitions,  the  Raman  vibration-rotation  bands,  which  are  displaced  in 
magnitude  from  the  laser  source  wave  number  by  >  1000  cm1,  can  be  observed.  Unlike 
pure  rotational  structures,  these  bands  provide  a  high  degree  of  molecular  specificity. 
Since  the  intensity  of  each  band  is  directly  proportional  to  the  molecular  species  number 
density,  the  intensity  of  the  N2  vibration-rotation  band  Q-branch  at  8283  A  can  be 
expressed  as 


1 A  Q)  = 


n(N2) 


CF(N2, 


Q)C,(T)l-] 


(6) 


Again,  if  temperature  has  been  measured  and  if  the  rotational  and  vibrational  modes  are 
in  equilibrium,  then  the  N2  number  density  can  be  determined  with  Eq.  (6).  The  plot  of 
Cf4(T)  as  a  function  of  temperature  is  given  in  Fig.  16. 

3.2  CALIBRATIONS 

Calibrations  were  conducted  before  each  test  period.  A  slow,  steady  flow  of  either 
pure  N2  or  air  maintained  atmospheric  pressure  in  the  test  cell.  These  calibrations  were 
used  to  determine  Cp(N2)  and  CF(air).  The  relative  sensitivity  of  the  detectors  at  7010  A 
(Cf  r  )  was  also  determined. 

Because  of  the  need  to  perform  radial  surveys  of  the  flow  field,  it  was  also  necessary 
to  calibrate  the  optical  sensitivity  of  the  system  as  a  function  of  the  radial  position  of 
the  laser  focal  volume.  With  a  flow  of  either  pure  N2  or  air,  the  Raman  signal  at  7010  A 
was  measured  over  the  lower  half  of  the  flow  field , at  0.3-cm  (1/8-in.)  intervals.  Only 
minor  changes  in  relative  sensitivity  were  noted  until  within  0.6  cm  (1/4  in.)  of  the  cell 
wall,  at  which  point  collection  solid  angle  began  to  be  rapidly  decreased.  Therefore,  no 
test  measurements  were  made  closer  to  the  wall  than  0.6  cm. 
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3.3  DATA  REDUCTION 

Nominally  fifteen  laser  pulses  were  used  for  each  measurement  point  within  the  flow 
field.  For  each  laser  pulse  an  intensity  ratio,  Rj,  was  determined  with  the  relation 


R*  = 


lmi(7061.2  A)  CFH 


I  .(7010  A) 

mi 


(7) 


and  a  corresponding  temperature,  Tj,  was  then  determined  by  graph.  The  ensemble 
average  of  these  temperature  values  provided  an  initial  value  of  average  temperature  at 
the  measurement  point  using 

_  k  x 

Tr  <R>  =  £  ~  (8) 

i-1  k 


where  k  is  the  number  of  laser  pulses  per  data  point.  The  fractional  rms  deviation  of  the 
initial  average  temperature  was  determined  from 

[  I  arp*]  ‘/! 

s-  -  L-- _ J  C9) 

T  k  Tr(R) 

in  which 


5Ti  =  T.  _  Tn(R) 

Because  of  the  nonlinearity  of  the  R  versus  T  curves  and  the  limited  number  of  laser 
pulses  per  data  point,  an  additional  correction  procedure  must  be  adopted  to  arrive  at  a 
true  average  temperature  value,  TR  (Refs.  8  and  9).  The  true  average  temperature  can  be 
estimated  as 


in  which 


(10) 


(11) 


13 


A  EDC-TR-79-74 


The  Mj  values  are  fluctuation  measures  defined  as  follows: 


M.  =  1  +  M' 

1  i 


(12) 


in  which  Mj  is  a  factor  to  account  for  detector  shot  noise  and  turbulent  medium 
fluctuations.  Mj  is  a  random  number  which  is  Poisson  distributed  about  unity;  i.e., 

M.  -  1,  S'  -  0  (13) 


To  evaluate  eM  it  must  be  noted  that  the  shot  noise  and  turbulent  medium  fluctuations 
are  uncorrelated,  a  condition  that  allows  the  neglect  of  the  product  term  of  M,  values. 
Further, 

4  -  (■“;/  +  K?  U4) 

in  which  Sr  is  the  fractional  rms  fluctuation  of  the  measured  intensity  ratio  value.  Sr  is 
calculated  from  the  following  equation: 

k  o 

y  (SRj)2 

i  =  1  k 

oR  =  - - - 

fl 

in  which 


and 


R 


SR.  =  Rj-R 


The  value  of  R  corresponding  to  Tr(R)  is  defined  to  be  R0.  For  300  <  Tr(R) 

<1200  K, 

d2T 

— -  -  -0.825554  x  105  +  (1.362582  x  106)  R  (16) 

dR2 


-(7.289328  x  106)  R2  +  (12.50538  x  106)  R3 
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For  1200  <  Tr(R)  <3000  K, 


rl^  X  ■! 

— -  =  1.152546  x  103  -(0.610152  x  104)  R 
dH2 


+  (3.473112  x  in4)  R2  -  (2.20656  x  104)  R3 


07) 


The  second-order  derivatives  were  obtained  from  a  Fifth-order  polynomial  fit  to  the  TR 
versus  R  curve  in  Fig.  12. 

Initial  air  number  densities  were  obtained  from 


_ 0> 

n(air) 


Cp(air)  Cp  (T  p)  O{7010  A) 


(18) 


However,  as  can  be  seen  in  Fig.  15,  CfjCD  is  not  a  linear  function  of  temperature. 
Therefore,  a  correction  is  required  to  obtain  the  true  number  density  values: 

_  _ (I) 

n(air)  =  n(air)  —  (n  (19) 

in  which 

( S T) 2  Cp(air)  1^  (7010  A)  (20) 

and  ' 

ST  =  ST(k)‘^  (21) 


d2C, 


dT- 


A  radial  profile  of  H2  number  density  was  obtained  at  one  axial  position.  However, 
because  of  a  lack  of  test  time  for  obtaining  a  calibration  in  H2,  the  density  values  had  to 
be  presented  as  only  relative  density  values.  At  the  same  axial  position  a  radial  profile  of 
N2  number  density  was  obtained  by  using  the  N2  vibration-rotation  band  and  the 
following  relation 


(22) 
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No  correction  was  applied  because  of  the  near-linear  dependence  of  Cf4  on  temperature 
over  the  range  of  temperatures  and  temperature  fluctuations  encountered. 


4.0  DISCUSSION 

4.1  PRESENTATION  OF  DATA 

Tabulations  of  the  laser-Raman  results  are  given  in  Tables  1  and  2.  Figures  17 
through  23  show  the  radial  variations  of  the  temperature  and  number  density  values; 
Figs.  24  and  25  show  the  axial  variations  of  number  density. 

Average  minimum  temperatures  are  approximately  320  K,  and  average  maximum  to 
temperatures  are  approximately  1500  K.  It  was  anticipated  that  maximum  temperatures 
near  2500  to  3000  K  would  be  observed  in  the  combustion  regions  of  the  flow,  and 
instantaneous  values  of  temperature  were  observed  to  be  this  high  (Figs.  26  and  27).  It  is 
speculated  that  the  high-temperature  combustion  zones  may  have  had  a  spatial  extent 
smaller  than  the  cylindrical  Raman  observation  volume  of  approximately  1-mm  diam  and 
3.1 8-mm  length. 

Number  density  values  were  obtained  predominantly  from  the  intensity 
measurements  of  the  rotational  structure.  The  N2/02  mole  fraction  ratio  was  unknown  a 
priori;  therefore,  number  density  values  were  reduced  in  accordance  with  the  assumption 
of  an  N2/02  mole  fraction  ratio  of  3.73  (i.e. ,  air).  This  assumption  results  in  a  large 
uncertainty  (discussed  in  Section  4.2)  for  the  combustion  zones  in  which  the  02 
concentration  is  depleted.  However,  for  regions  at  or  near  the  axial  centerline  and  near 
the  nozzle  exit  plane,  the  assumption  is  valid,  and  it  should  be  noted  that,  as  expected, 
the  number  densities  are  near  atmospheric  values. 

4.2  DATA  UNCERTAINTY 

Data  uncertainty  values  are  given  for  each  data  point  in  Tables  1  and  2.  These  were 
calculated  with  the  relation 

U  =  B  +  S  (23) 

in  which  B  represents  the  fractional  systematic  error  and  S  represents  the  fractional 
random  error.  Equation  (23)  gives  a  range  of  uncertainty  (minimum,  maximum)  for  the 
measurements.  For  the  temperature  measurements,  the  uncertainty  was  calculated  with 

Urp  =  By  +  Sy  (24) 
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The  values  are  also  given  for  each  data  point  in  Tables  1  and  2.  The  systematic  error 
(or  bias)  results  from  the  assumption  that  the  N2/O2  mole  fraction  was  3.73  (i.e.,  the 
same  as  for  air).  As  shown  in  Fig.  1 2,  the  maximum  uncertainty  would  result  if  the  O2 
were  completely  depleted  and  if  the  R  versus  TR  curve  for  N2  had  been  used.  By  using 
the  results  shown  in  Fig.  1 2,  the  maximum  systematic  error  was  calculated  as  a  function 
of  temperature  and  is  plotted  in  Fig.  8.  From  this  plot,  the  values  of  BT  were 
determined. 

The  uncertainty  for  the  density  measurements  was  calculated  using 

-  -  B„  i  s-  (25) 

The  systematic  error  results  from  the  N2/O2  mole  fraction  assumption.  As  shown  in  Fig. 
15,  the  maximum  uncertainty  would  result  if  the  02  were  completely  depleted  and  if  the 
CF2(T2)  versus  TR  curve  had  been  used.  The  maximum  systematic  error  calculated  with 
the  use  of  Fig.  1  5  is  plotted  as  a  function  of  temperature  in  Fig.  28.  From  this  plot  the 
values  of  Bn  were  determined. 

The  random  error  of  the  density  measurements  was  calculated  as  follows: 

S5  =  the  fractional  rms  deviation  of  the  intensity  measurement 

£  (aF  .(7010A) 

s-  =  L=i^! - 

n  - 5 — 

(k)  lc  (7,010  A) 

m 

SJf  =  the  fractional  deviation  of  the  CFl(TR)  value  resulting  from  the  imprecision  of 
the  temperature  measurement 


V, 

(26) 


Si7  =  ST  TIi 


(27) 


For  300  K  <  TR  <  1 200  K, 


— 1 -  =  -2.803  029225  x  10“  2  +  (1.357945891  x  10“  4)  T 

dT 

-  (2.468528696  x  10  “7)  T2  +  (1.975217603  x  10“  10)  T3 

-  (5. 840266985  x  10 “  14)  T4  (28) 
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For  1200  K  <  TR  <3000  K, 

<K^>  =  7.849131548  x  10“4  -  (1-61032032  x  10“6)  T 

dT 

+  (1.218008693  x  10“9)  T2  -  (3.981211523  x  10“ 13)  T3 
+  (4.783911592  x  10“ 17)  T4  '  (29) 

S„  =  the  random  error  of  the  density  measurements  as  follows: 

S_.  [<sj>2  +  (Sf).2]  %  (30) 


4.3  COMPARISON  OF  RESULTS  FOR  THE  TWO  TEST  CASES 

The  contrast  between  the  flow  fields  of  the  two  test  configurations  is  striking.  For 
DTINDR  =  2,  the  fuel-to-air  mass  ratio  (F/A)  was  0.00345,  which  is  quite  lean,  and 
apparently  the  primary  air  flow  dominates  the  flow  structure.  As  observed  in  Figs.  17 
and  18,  the  combustion  is  limited  to  the  radial  positions  greater  than  1  in.  and  to  axial  \ 
regions  within  10.5  in.  of  the  nozzle  exit  plane.  The  inner  air  core  remains  relatively  cool 
until  axial  positions  beyond  21  in.  are  reached,  at  which  point  it  appears  that  the  hot 
combustion  gases  and  the  air  flow  mix  together  to  achieve  a  uniform  gas  temperature  of 
—550  K.  For  DTINDR  =  10.  the  F/A  was  increased  to  0.107.  The  combustion  process 
begins  quite  early  for  this  configuration  (axial  position  <  2.5  in.)  and  very  close  to  the 
boundary  between  the  primary  air  flow  and  the  H2  flow.  The  primary  air  core  remains 
cool  until  an  axial  position  of  5.25  in.;  however,  beyond  this  position,  the  temperature  of 
the  central  core  progressively  rises  as  the  combustion  process  continues  and  shifts  to  the 
central  region  of  the  flow.  In  contrast  to  the  DTINDR  =  2  configuration,  average 
temperatures  on  the  centerline  rise  to  nearly  1400  K. 

The  axial  centerline  variation  of  air  number  density  (see  Fig.  24)  also  indicates  a 
significant  difference  between  the  two  test  configurations.  For  DTINDR  =  2.0,  the  air 
number  density  on  the  axial  centerline  remains  nearly  constant  out  to  the  24-in.  axial 
position,  at  which  point  it  begins  to  decrease  as  centerline  temperatures  begin  to  increase. 

In  contrast,  the  axial  centerline  air  number  density  for  DTINDR  =  10.0  rapidly  decreases 
with  axial  distance  from  the  nozzle  exit. 
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4.4  COMPARISON  OF  THE  RESULTS  WITH  CALCULATIONS 

Equilibrium  chemistry  calculations  were  available  only  for  the  DTINDR  =  10.0 
configuration.  As  shown  in  Figs.  19  and  20,  the  measured  temperature  values  disagree 
markedly  with  the  equilibrium  chemistry  predictions.  At  the  axial  position  of  2.5  in.,  a 
peak  average  temperature  of  1900  K  is  predicted  within  0.25  in.  of  the  axial  centerline. 
Measurements  show  a  peak  temperature  of  1100  K  at  the  0.6-in.  radial  position.  Similar 
disagreements  are  observed  in  other  axial  positions.  However,  the  measured  temperature 
profile  at  the  26.25-in.  axial  position  indicates  that  measurements  and  predictions  for 
large  axial  positions  are  beginning  to  move  toward  agreement.  Furthermore,  the  measured 
air  number  densities  at  the  21-  and  26.25-in.  axial  positions  bracket  the  predictions  of  N2 
number  density  as  shown  in  Fig.  23. 


5.0  SUMMARY 

The  laser-Raman  technique  has  been  successfully  used  to  partially  map  the  flow  field 
of  an  H2/air  combustor  for  two  widely  different  test  conditions.  The  system  performed 
admirably  under  what  must  truly  be  described  as  hostile  conditions.  The  test  cell  area 
was  hot  and  humid  (as  a  result  of  the  water  spray  coolers)  during  operation  as  well  as 
dirty  and  dusty  at  all  times.  As  a  result  of  energy  conservation  measures,  the  equipment 
was  frequently  exposed  to  subfreezing  temperatures  when  not  in  use.  Furthermore, 
during  one  test  period,  the  viewport  was  blown  out  twice  during  combustion  ignition,  but 
the  optical  system  was  only  slightly  damaged.  As  there  were  many  possibilities  for 
equipment  malfunction  to  result  directly  from  the  test  cell  area  environment,  it  is  a 
credit  to  the  durability  of  the  system  that  only  moderate,  easily  corrected, 
electro-optical-mechanical  difficulties  were  encountered. 

The  original  experiment  plan  included  mapping  the  flow-field  H2  density  and  N2 
density  using  the  H2  rotational  lines  and  the  N2  vibration-rotation  band,  respectively.  In 
addition  the  H2  rotational  lines  were  to  be  used  to  provide  more  accurate  determination 
of  temperature  in  regions  of  high  H2/low  N2  concentrations.  However,  limited  resources 
prevented  making  these  measurements.  As  a  result,  the  measurement  uncertainties  are 
much  higher  than  expected,  and  this  limits  the  usefulness  of  the  data  for  combustion 
modeling.  However,  a  further  analysis  of  the  data  [i'.e.,  comparisons  with  flow-field 
parameter  calculations,  pressure  probe  measurements,  infrared  (IR)  measurements,  and 
laser  velocimeter  measurements!  can  reduce  the  uncertainty  of  the  Raman  measurements 
by  providing  a  more  accurate  estimate  of  the  N2/02  mole  traction  ratio.  Although  such 
analysis  is  beyond  the  scope  of  this  report,  the  following  conclusions/observations  can  be 
made  for  these  experiments: 
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1.  One  of  the  main  impediments  to  faster  data  acquisition  was  the  poor 
performance  of  the  hydraulic  scanning  table  on  which  the  laser-Raman 
optical  system  was  mounted.  This  table  should  be  replaced  by  one  that  is 
sturdier  and  more  capable  of  easy,  accurate,  and  rapid  positioning. 

2.  Background  radiation  from  the  flow  field  was  undetectable  because  of  the 
cleanliness  of  the  flow  Field.  The  absence  of  particulates  also  eliminated 
laser-induced  particulate  incandescence,  which  would  have  been  manifested 
by  off-band  radiation  during  laser  firings. 

3.  In  the  high-temperature  zones  of  the  flow  fields,  instantaneous 
temperature  values  approached  the  expected  peak  level  of  —2500  K. 
However,  average  temperature  values  were  never  higher  than  1500  K.  The 
spatial  resolution  of  the  measurements  may  not  have  been  sufficient  to 
resolve  the  peak  temperature  zones.  Furthermore,  if  the  turbulence 
frequency  reaches  the  order  of  kilohertz,  then  the  number  of  samples  used 
for  a  data  point  must  be  increased  by  nearly  two  orders  of  magnitude  to 
obtain  a  more  reliable  average  temperature. 

4.  For  the  DTINDR  =  2  configuration,  the  peak  combustion  zone  is  near  the 
5.21-in.  position,  and  the  highest  average  temperature  (~1400  K)  occurs 
there  near  the  test  cell  wall.  Axial  centerline  temperatures  begin  to 
increase  only  at  the  21.07-in.  axial  position,  and  the  flow  appears  to  be 
uniform  with  a  nominal  temperature  of  550  K  by  the  31.5-in.  axial 
position. 

5.  The  peak  combustion  zone  is  much  more  difficult  to  define  for  the 
DTINDR  =  10  configuration.  Temperatures  as  high  as  1100  K  were  found 
at  the  2.59-in.  axial  position  and  within  1  in.  of  the  axial  centerline 
followed  by  decreasing  temperatures  toward  the  test  cell  wall.  Centerline 
temperatures  rapidly  increased  as  axial  distance  from  the  nozzle  increased; 
even  at  the  26.27-in.  axial  position,  the  flow  core  was  at  1400  K,  with 
temperatures  decreasing  toward  the  test  cell  wall. 

6.  For  the  DTINDR  =  2  configuration,  the  number  density  radial  distribution 
at  the  0.5 5-in.  axial  position  shows  the  broad  uniform  core  of  the  air  flow 
corresponding  to  the  1.31 -in.  radius  of  the  primary  nozzle.  Beyond  the 
1.31-in.  radial  position,  the  number  density  profile  exhibits  a  “sawtooth” 
behavior  that  corresponds  to  the  secondary  flow  hole-separation  in  the 
baffle  plate. 
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7.  For  the  DTINDR  =  10  configuration,  the  nearest  axial  position  at  which  a 
radial  scan  was  made  was  2.6  in.  At  this  position  the  H2  secondary  flow 
intruded  into  the  primary  air  flow.  Even  so,  the  primary  air  flow'  core 
corresponding  to  the  0.262-in.  radius  of  the  primary  nozzle  is  w ell-defined. 
Within  0.5  in.  of  the  axial  centerline,  the  H2  concentration  increased  by 
nearly  two  orders  of  magnitude. 

8.  The  general  trend  of  the  number  density  radial  profiles  for  the  DTINDR  = 
2  configuration  is  a  broad,  uniform  air  core  with  decreasing  number 
density  toward  the  test  cell  wall.  For  the  DTINDR  =  10  configuration,  the 
number  density  scans  at  axial  positions  of  10.44  in.  or  less  have  a 
U-shaped  profile.  That  is,  the  number  densities  decrease  from  the  central 
core  value,  bottom  out,  and  then  gradually  increase  toward  the  test  cell 
wall. 

9.  Equilibrium  chemistry  calculations  for  the  DTINDR  =  10  configuration 
fail  to  predict  the  observed  temperature  and  density  profiles  for  axial 
positions  of  15.75  in.  or  less.  However,  for  the  21.02-  and  26.27-in.  axial 
positions,  the  measurements  and  predictions  agree  to  within  ±20  percent. 
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Up  is  the  primary  jet  velocity  at  nozzle  exit. 
ug  is  the  secondary  jet  velocity  at  nozzle  exit. 

Mean  streamline  pattern  shown  in  upper  half  figure. 

Mean  axial  velocity  profiles  shown  in  lower  half  figure. 

Figure  1.  Turbulent,  ducted,  mixing  system  with  recirculation. 
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Figure  2.  Recirculating  flow  combustor  simulation  facility 
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a.  Data  acquisition  electronics 

Figure  4.  Block  diagram  of  electronic  systems  applicable  during  the  data 
acquisition  interval,  Q-switch/integrator  technique. 
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b.  Laser  system  electronics 
Figure  4.  Concluded. 
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Figure  6.  Pure  rotational  raman  spectrum,  N2  and  02  at  T  =  300  K; 
[N2  ]/[02  ]  =  3.73. 


A  E  DC-TR-79-7 4 


Intensity 


6900  6950  7000  7050  7100  7150  7200 

Wavelength,  S 


Figure  7.  Pure  rotational  Raman  spectrum,  l\l2  and  02  at  T  =  900  K; 
[^2  ]/[02  ]  =  3.73, 


A  EDC-TR  -79-74 


Intensity 


Wavelength , 


Figure  8.  Pure  rotational  Raman  spectrum,  N2  and  02  at  T  =  1500  K; 
[N2  I/[02  ]  =  3.73. 
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Figure  9.  Pure  rotational  Raman 
[N2  ]/[H2  ]  =  0.667. 
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Figure  10.  Pure  rotational  Raman  spectrum,  N2  and  H2  at  T  =  900  K; 
[N2}/[H2]  =  0.667. 
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Figure  11.  Pure  rotational  Raman  spectrum,  l\J2  and  H2  at  T  =  1500  K; 
[N2]/[H2]  =  0.667. 
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Figure  17.  Isometric  plot  of  measured  radial  temperature  distribution, 
DTINDR  =  2;  0.55-,  5.25-,  10.5-in.  axial  positions. 
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Figure  19.  Isometric  plot  of  measured  radial  temperature  distribution 
DTINDR  =  10;  2,5-,  5.25-,  10.5-in.  axial  positions. 
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Figure  20.  Isometric  plot  of  measured  radial  temperature  distribution, 
DTINDR  =  10;  15.75-,  21-,  26.25-in.  axial  positions. 
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Figure  21.  Air  number  density  as  a  function  of  flow-field  position,  DTINDR  =  2. 
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Figure  22.  Number  density  as  a  function  of  flow-field  position, 
DTINDR  =  2. 
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Figure  23.  Number  density  as  a  function  of  flow-field  position, 
DTINDR  =  10. 
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Figure  25.  Axial  variation  of  n (air)  for  both  test  configurations. 
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Figure  26.  Instantaneous  temperatures  versus  laser  pulse  number. 
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Figure  27.  Instantaneous  temperatures  versus  laser  pulse  number. 
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Table  1.  Measured  Values  of  Temperature  and  Number  Density  and  Estimated 
Uncertainty,  DTINDR  =  2 
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Table  1.  Concluded 
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Table  2.  Measured  Values  of  Temperature  and  Number  Density  and  Estimated 
Uncertainty,  DTINDR  =  10 


Test  Period/ 
Data  Point 
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Position, 
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2.0  x  10 


1.3  x  10 


1.2  x  10 


1.4  x  10 


-0.33,  -0.26 
-0.44,  -0.26 
-0.226,  -0.468 
-0.247,  -0.379 
-0.212,  -0.420 
-0.268,  -0.408 
-0.308,  -0.415 
-0.297,  -0.387 
-0.305,  -0.385 
-0.291,  -0.393 
-0.206,  -0.340 
-0.279,  -0.431 
-0.298,  -0.510 
-0.371,  -0.501 
-0.361,  -0.451 
-0.237,  -0.535 
-0.357,  -0.433 
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Table  2.  Continued 


Test  Period/ 
Data  Point 

Axial 

Position, 

in. 

Radial 

Position, 

in. 

Tr(R). 

K 

*1 

V 

K 

UT 

n(air) , 
on'3 

S5 

Un 

6/19 

21 

02 

1.50 

1203 

0.040 

1197 

0.090,  0.170 

8.0  x  1017 

0.092 

-0.290,  -0.474 

6/20 

1.74 

1140 

0.062 

1073 

0.065,  0.189 

1.3  x  1018 

0.038 

-0.316,  -0.392 

6/21 

1.99 

976 

0.042 

981 

0.078,  0.162 

8.7  x  1017 

0.062 

-0.268,  -0.342 

6/22 

2.26 

1078 

0.037 

1073 

0.090,  0.164 

1.3  x  1018 

0.023 

-0.329,  -0.375 

6/23 

2.50 

738 

0.053 

740 

0.051,  0.157 

2.5  x  1018 

0.048 

-0.207,  -0.303 

6/24 

26 

27 

0.0 

1362 

0.044 

1351 

0.096,  0.184 

1.2  x  1018 

0.061 

-0.353,  -0.475 

6/25 

0.26 

1421 

0.041 

1420 

0.099,  0.181 

1.8  x  1018 

0.070 

-0.356,  -0.496 

6/26 

0.49 

1335 

0.064 

1315 

0.076,  0.204 

7.7  x  1017 

0.112 

-0.294,  -0.519 

6/27 

0.99 

1403 

0.047 

1390 

0.089,  0.191 

1.8  x  1018 

0.040 

-0.382,  -0.462 

6/28 

1.51 

1142 

0.051 

1094 

0.078,  0.180 

1.3  x  1018 

0.030 

-0.328,  -0.388 

6/29 

2.0 

1060 

0.055 

1058 

0.072,  0.182 

9.1  x  1017 

0.047 

-0.303,  -0.397 

6/30 

2.5 

809 

0.039 

809 

0.072,  0.150 

3.3  x  1018 

0.058 

-0.219,  -0.335 

6/31 

2 

59 

0.0 

397 

0.040 

401 

0.030,  0.110 

1.5  x  1019 

0.040 

-0.052,  -0.132 

6/32 

5 

29 

0.0 

446 

0.032 

449 

0.043,  0.107 

1.0  x  1019 

0.044 

-0.079,  -0.167 

6/33 

10 

49 

0.0 

599 

0.070 

598 

0.039,  0.179 

5.6  x  1018 

0.049 

-0.151,  -0.249 

6/34 

0.48 

901 

0.062 

907 

0.058,  0.182 

4.2  x  1018 

0.078 

-0.227,  -0.383 

6/35 

1.01 

897 

0.042 

900 

0.078,  0.162 

2.8  x  1018 

0.057 

-0.247,  -0.361 

6/26 

_ 

1.51 

980 

0.039 

986 

0.081,  0.159 

2.3  x  1018 

0.037 

-0.294,  -0.368 
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Table  2.  Continued 


Axial 

Position, 

in. 

Radial 

Position, 

in. 

*T 

V 

K 

UT 

n(air), 

-3 

cm 

s- 

n 

U 

n 

Test  Period/ 
Data  Point 

Tr(R), 

K 

5/3 

2.1 

>5 

0.0 

357 

0.040 

360 

0.000,  0.080 

2.0 

X 

1019 

0.044 

-0.018,  -0.106 

5/4 

0.13 

310 

0.036 

312 

-0.030,  0.042 

2.2 

X 

1019 

0.056 

0.051,  -0.061 

5/5 

0.26 

570 

0.070 

589 

0.028,  0.168 

4.6 

X 

io18 

0.039 

-0.153,  -0.236 

5/6 

0.39 

657 

0.136 

698 

-0.036,  0.236 

2.5 

X 

io18 

0.067 

0.039,  -0.201 

5/7 

2.63 

0.49 

810 

0.096 

810 

0.015,  0.207 

2.1 

X 

1018 

0.066 

-0.212,  -0.344 

5/8 

0.63 

1070 

0.134 

1044 

-0.010,  0.258 

2.2 

X 

io18 

0.081 

-0.265,  -0.427 

5/9 

0.74 

1106 

0.088 

1050 

0.036,  0.212 

2.2 

X 

io18 

0.055 

-0.292,  -0.402 

5/10 

0.88 

1060 

0.049 

1055 

0.076,  0.174 

2.1 

X 

io18 

0.071 

-0.278,  -0.420 

5/11 

1.01 

1021 

0.087 

1052 

0.037,  0.211 

1.9 

X 

io'8 

0.072 

-0.276,  -0.420 

5/12 

1.11 

955 

0. 108 

1009 

0.013,  0.229 

1.7 

X 

io18 

0.064 

-0.273,  -0.401 

5/13 

1.24 

663 

0.  104 

706 

-0.004,  0.204 

2.5 

X 

io18 

0.054 

-0.191,  -0.299 

5/14 

1.37 

763 

0.  138 

769 

-0.033,  0.243 

2.3 

X 

io18 

0.  105 

-0.162,  -0.372 

5/15 

1.61 

598 

0.  115 

686 

-0.015,  0.215 

2.1 

X 

io'8 

0.067 

-0.170,  -0.304 

5/16 

1.84 

570 

0.  113 

622 

-0.015,  0.215 

2.5 

X 

io18 

0.060 

-0.150,  -0.270 

5/17 

2.13 

459 

0.063 

486 

0.015,  0.141 

3.6 

X 

io’8 

0.052 

-0.094,  -0.198 

5/18 

2.62 

2.35 

488 

0.061 

503 

0.019,  0.141 

7.2 

X 

io18- 

0.038 

-0.112,  -0.188 

5/19 

2.62 

0.0 

361 

0.039 

361 

0.001,  0.079 

2.1 

X 

io19 

0.015 

-0.045,  -0.075 

5/20 

5.25 

0.0 

375 

0.068 

387 

-0.007,  0.129 

1.3 

X 

1019 

0.047 

-0.036,  -0.130 

5/21 

0.14 

395 

0.060 

405 

0.010,  0.130 

1.1 

X 

io’9 

0.076 

-0.017,  -0.169 

5/23 

0.38 

801 

0.133 

800 

-0.023,  0.243 

3.8 

X 

io18 

-0.225,  -0.325 

5/24 

0.50 

930 

0.097 

946 

0.023,  0.217 

3.6 

X 

io18 

0.067 

-0.193,  -0.327 

5/25 

0.62 

810 

0.107 

810 

0.004,  0.218 

3.1 

X 

io18 

0.074 

-0.206,  -0.354 
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Table  2 


Test  Period/ 
Data  Point 

Axial 

Position, 

in. 

Radial 

Position, 

in. 

Tr(r), 

K 

s? 

V 

K 

mm 

5.25 

0.75 

904 

0.079 

914 

0.86 

950 

0.067 

960 

5/28 

1.11 

840 

0.080 

841 

5/29 

1.38 

689 

0.083 

709 

5/30 

5/31 

1.63 

773 

0.067 

843 

1.87 

693 

0.087 

715 

5/32 

2.13 

686 

0.032 

709 

5/33 

2.38 

640 

0.027 

662 

5/34 

’ 

0.0 

331 

0.038 

334 

5/36 

10.44 

0.01 

480 

0.040 

488 

5/37 

0.11 

520 

0.064 

552 

5/38 

0.25 

590 

0.107 

661 

5/39 

0.38 

630 

0.054 

649 

5/40 

0.50 

640 

0.049 

654 

5/41 

0.64 

840 

0.048 

841 

5/42 

0.77 

940 

0.053 

951 

5/43 

10 

43 

0.88 

850 

0.069 

853 

5/44 

1.12 

920 

0.044 

927 

5/45 

1.35 

870 

0.043 

872 

5/46 

1.64 

750 

0.032 

750 

5/47 

1.88 

830 

0.045 

830 

5/48 

2.14 

840 

0.033 

840 

5/49 

1 

1 

0.01 

450 

0.054 

480 
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NOMENCLATURE 

B  .  Fractional  systematic  error 

Bn>  Bt  Fractional  systematic  error  of  number  density  and  temper¬ 

ature  measurements,  respectively 


Cp(air),  CF(H2},  Cj- C^2 )  Calibration  factors  for  air,  H2,  and  N2  density  determination 

using  rotational  line  intensities 

Cp(N2l  Q)  Calibration  factor  for  N2  density  determination  using  the  N2 

■  vibration-rotation  band  Q-branch 


cFl(T),  c,,2(T) 


Relative  temperature  dependence  of  the  intensity  of  the 
rotational  lines  at  7010  A  for  air  and  pure  N2,  respectively 


C,.3(T) 


Relative  temperature  dependence  of  the  intensity  of  the  J=l, 
H2  rotational  line 


Cf-  4(T) 


Relative  temperature  dependence  of  the  intensity  of  the  N2 
vibration-rotation  band  O-branch 


Q.  r 


Relative  sensitivity  of.  the  two  detection  channels 


DT1NDR 


Duct-to-inner-nozzle  diameter  ratio 


F/A 


Fuel-to-air  mass  ratio 


1(7010.  A),  1(7061.2  A) 


I(J=1,  H2 ),  I(J=2,  H2) 


Im 

i(n2,  9) 

j 

k 


Intensity  of  the  rotational  Raman  lines  at  7,010  A  and 
7061.2  A,  respectively 

Intensity  of  the  J=1  and  J=2  H2  rotational  Raman  lines, 
respectively 

Intensity,  measured 

Intensity  of  the  N2  vibration-rotation  band  Q-band 
Intensity,  test  conditions,  measured 

Rotational  level  quantum  number 
Number  of  laser  pulses  used  per  data  point 
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M,  M' 
n(air) 
n(N  2 ) 
n(Ha) 

Q 

R,  R 

R(Tr) 

Ro 

S 

SH 


Sr 


St 


St 

T,  Tr 
Tr(R) 

Tr 


Fluctuation  measures 
Number  density  of  air 
Number  density  of  N2 
Number  density  of  H2 

Q-b  ranch 


Measured  intensity  ratio  and  average  measured  intensity  ratio, 
respectively 

Calculated  rotational  line  intensity  ratio 
Intensity  ratio  corresponding  to  Tr(R) 


Fractional  random  error 

Fractional  random  error  of  the  average  number  density 
measurement 

Fractional  rms  deviation  of  intensity  values  used  to  determine 
average  number  density 

Fractional  deviation  of  CFl(TR)  values  resulting  from 
imprecision  of  the  Tr  determination  '  • 

Fractional  rms  deviation  of  the  intensity  ratios  used  to 
determine  an  average  temperature 

Fractional  random  error  of  the  average  temperature  measure¬ 
ment 

j 

Fractional  rms  deviation  of  the  temperature  values  used  to 
determine  an  average  temperature 

Temperature  and  rotational  temperature,  respectively 

Initial  average  temperature  value 

Final,  corrected  average  temperature  value 
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U;  UT ;  UN 


X 

51(7010  A) 

6R 

6T 


cm 


fn 

SUBSCRIPTS 

a,b 

c 

i 

m 

SUPERCR1PTS 


t 


Total  uncertainty;  total  uncertainty  in  temperature  measure¬ 
ment;  total  uncertainty  in  number  density  measurement 

Mole  fraction  of  species  X 

Measured  intensity  residuals 

Measured  intensity  ratio  residuals 

Measured  temperature  residuals 

Correction  term  for  obtaining  final  average  temperature 
Correction  term  for  obtaining  final  average  number  density 

Denotes  the  two  detector  channels 
Indicates  a  calculated  value 
indicates  an  instantaneous  value  i 

Indicates  a  measured  value 

Indicates  an  ensemble  average  of  the  quantity  beneath  the  bar 
Indicates  the  value  of  a  parameter  at  atmospheric  conditions 
Indicates  the  value  of  a  parameter  at  test  conditions 
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